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c/3 ' Sg is the most stable compound of elemental sulfur in solid and liquid phases, at 



ambient pressure and below 400K. Three crystalline phases of Sg have been clearly 



F^ ' identified in this range of thermodynamic parameters, although no calculation of its 

a' 

O ' phase diagram has been performed yet. a- and 7-S8 are orientationally ordered crys- 

tals while /J-Sg is measured as orientationally disordered. In this paper we analyze 

m 

K^ ' the phase diagram of Sg crystals, as given by a simple and fiexible molecule model, 

Os ■ 

iy~, ' via a series of molecular dynamics (MD) simulations. The calculations are performed 

o' 

in the constant pressure- constant temperature ensemble, using an algorithm that is 

5^ ' able to reproduce structural phase transitions. 

-(— > ' 
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'^ ' Introduction: 

^ ' Elemental sulfur shows a rich variety of molecular forms and therefore of chemical and 

physical properties in the crystalline and liquid phasestl'cl. Recently, a high pressure super- 

X' 

Vh ' conductor transition has been found, with Tc=17K (at 162GPa), the highest known for any 

elemental solidtl. Natural sulfur is of relevance in geophysics, astrophysics, material sciences 
and also of massive industrial uselll. Nevertheless, the complex phase diagram and properties 
of its condensed phases are far from well knownElo. 

The most stable sulfur allotrope at ambient temperature and pressure (STP) is Sg, a 
crown-shaped cyclic Sg molecule, stable in solid, liquid and gas phases. 



Orthorhombic a-Sg is the crystalline form stable at STP. Figs la shows its "crankshaft" 
structure with four molecules in the primitive unit cell. The space group is D^^ (Fddd), with 
16 molecules in this non-primitive orthorhombic cell. If a-Ss is slowly heated, it shows a 
phase transition to monoclinic pSs at 369K (Fig. lb), which melts at 393Ktl. Nevertheless, 
most a-Sg crystals do no easily convert to /^-Sg, they melt, instead, at 385.8KU'@. /5-Sg is a 
monoclinic crystal, with six molecules in the primitive unit cell, two of them orientationally 
disordered^. A third crystalline allotrope has been observed: 7-Sg (Fig. Ic), that can be 
obtained from solutions of Sg or from its melta'cl. The space group is C2/J (P2/c), with four 
molecules in the pseudo-hexagonal closed-packed unit celB. The density of this allotrope, 
at STP, is 5.8 % higher than that of a-Sg. 

In spite of its relevance, calculations on the crystalline phase diagram of elemental sulfur 
molecules are scarce. Previous molecular dynamics (MD) simulations of crystalline Sg were 
limited to constant-volume simulations at normal pressure and temperatures of a-SgQ, with 
the only aim to study the intramolecular frequencies obtained by Raman and ir. measure- 
ments. The same scope was followed by numerous works that measured crystalline frequen- 
cies and performed calculations based on lattice dynamics sums in different aproximationda. 
No attempt to reproduce the phase diagram of this molecule, including determination of 
crystalline structures, lattice dynamics, or thermodynamic properties of /?- and 7-Sg as a 
function of temperature was intended. Here we study the phase diagram of Sg crystals, as 
given by a simple model of a flexible molecule, by performing a series of constant pressure - 
constant temperature MD simulations, at several temperatures and zero pressure. 

The inter- and intramolecular potential model: 

The main problem in this type of calculation is the inter- and intramolecular potential 
model of this flexible molecule. The above mentioned constant volume MD simulations 
of a-Sg were performed with an intermolecular potential model of the atom-atom type, to 
study the splitting of intramolecular frequencies in a sample of 64 molecules of a -SgQ. In our 
constant temperature - constant pressure MD simulations, these atom-atom parameters are 



slightly changed to improve the fit of our results to the experimental unit cell volume and 
the estimated configurational energy of a-Sg at 300K and Okbar pressure. The parameters 
of our LJ potential model, for S-S interactions are: e= 1.70kJ/mol and cr=3.39 A. 

In a similar way, our intramolecular potential parameters for bending and torsional angles 
of the flexible Sg molecule are based on those of refH Our molecule model takes into account 
all intramolecular modes that mix with lattice modes and can therefore be relevant for 
structural phase transitions. As calculated in ref.Q and also measured in all crystalline phases 
at all temperatures, three main regions are clearly separated in the vibrational spectra: 
between 10 and ~100cm^^ are found the lattice and torsional modes, between 150 and 
250cm~^ the bending modes and from 400 to 500cm~^ the stretching modes, well above in 
energies than the rest of vibrational modes. The bond length is therefore held constant at 
2.060A in our calculations. Our flexible molecule model allows, instead, the mixing of the 
bending and torsional intramolecular motions with lattice modes in a straightforward way. 
The intramolecular potential for bending angles S-S-S is taken harmonic 

with a force constant of C/3=25725kB/rad^ and /5o= lOSdeg., k^ is the Boltzmann constant. 
The intramolecular potential for torsion (r) angles is a double well: 

V{t) = A^ + Br cos(r) + C^ cos^(r) + D^ cos^(r), 

with A^=57.192 k^, 5^=738.415 k^, ^=2297.880 k^ and D^=557.255 k^. These param- 
eters describe a double well with minima at r=l 98.8deg., and a barrier height of about 
9kJ/mol at r=180deg. The barrier height at r=Odeg. is of 30kJ/mol, out of the range of 
energies explored in these simulations. 

Calculations: 

The phase diagram and dynamical properties of Sg crystals, as given by this molecular 
model, are studied in the (N,P,T) ensemble, at zero pressure, by a series of classical MD 



simulations. The external pressure is considered isotropic and the MD algorithm allows 
volume and shape fluctuations of the sample, in order to balance the internal stresses with 
an external isotropic pressure^. The algorithm controls pressure via an extended system 
which includes as extra variables the MD box parameters. The temperature control of the 
sample follows the approach of Nosel2ry, which also includes an external variable. The 
equations of motion of these flexible molecules are integrated using the Verlet algorithm for 
the atomic displacements and the Shake algorithm for the constant bond length constraints 
on each moleculec^ilfl. The flnal MD algorithm is totally similar to that used in a study of 
black Newton fllmsEj. 

As the starting point of our simulations, we considered the available structural data at 
STP. The molecules were located at the experimental sites, with the measured orientation 
but with their isolated D4d symmetry. In the course of the simulation, the molecules distort 
to attain, on a time average, their site symmetry. Following this procedure, the only problem 
found was with the /J-Sg sample. In this case the calculations were initiated with the LJ 
parameter e multiplied by a factor that in a few (about 100) time steps was varied from 0.1 
to 1. This procedure allowed the molecule to attain its initial distortion in the volume of 
the experimental unit cell. 

In a flrst step we stabilized our samples, at the temperatures quoted in Table I, with a 
constant volume algorithm for 30000 time steps (each one of O.Olps.). Later on, the runs 
in the (N,P,T) ensemble were performed by lowering and increasing the temperature (up 
to 500K) in steps of 25 or 50K. At each point of the phase diagram the samples are sta- 
bilized for 20000 to 30000 time steps and measured in the following lOps. Near the phase 
transitions the stabilization times were increased several times. The sample of a-Sg crystals 
consisted of 3*3*2 orthorhombic cells, a sample of 288 molecules. Some points of the phase 
diagram were recalculated with a sample of 4*4*2 cells (512 molecules). The sample of jS-Ss 
crystals consisted of 4*4*4 cells (384 molecules) and that of 7-S8 crystals of 5*3*5 cells (300 
molecules) . 



Results: 

a) a-Sg crystals: 

a-Ss crystallizes, at STP, in the ortliorhombic space group F^ (D2h) with 16 molecules in 
a face-centered unit cell, the primitive cell contains 4 molecules at C2 sites. Fig. la includes 
part of our sample at 300K and shows the structure in the [110] direction. The molecular 
rings lie parallel to the crystallographic c axis and are alternately oriented parallel to (110) 
and (110) crystalline planes. Table I shows the experimental lattice parameters and our 
calculated values, the largest difference is found in the value of b axis, 2 % larger than its 
experimental value. The molecular array and their corresponding orientations follows the 
experimental data. Our calculated configurational energy is -103.5(2)kJ/mol at 300K, from 
which the heat of sublimation can be estimated as lOl.OkJ/mol and compares well with a 
measured value of 101.8(2)kJ/molt3. Also our calculated pair distribution function g2(r), 
not included here, compares well with neutron measurements of crystalline a-Sg at SOOKta. 

Fig. 2 shows, for a—, (3— and 7 — 5*8, the calculated configurational energy and volume 
per molecule as a function of temperature. Our simulations reproduce the experimental 
data of the three crystals at STP. Unfortunately, our a-Sg sample of 288 molecules turns 
out to be unstable for T<200K, and distorts to a monoclinic cell, with molecular locations 
and orientations closely resembling those of a-Ss- This result was later on confirmed with 
a larger sample of 512 molecules (4*4*2 orthorhombic cells). Below 200K this large sample 
does not distort to monoclinic, but nevertheless its configurational energy is higher, and its 
volume per molecule is larger, than those calculated for the sample of 288 molecules. Also, 
the fluctuations of the energy and volume, of the sample of 512 molecules, are larger for 
T<200 than their value at 200K. 

The coincidence of the calculated distortion at 200K and the experimental fact that 
crystals of pure Sg are greenish-yellow at 293K and turn white below 193K, made us review 
the experimental data on this phase. Usually the reference quoted for lattice parameters as 
a function of temperature is ref.li3. These values were obtained using X-ray powder diffrac- 
trometry, on different samples and using only three lines. Our calculated powder diffraction 



patterns, not included here, were similar for the orthorhombic and monoclinic samples and 
the scarce experimental data could be not enough to detect the changes. Nevertheless, there 
is a careful single crystal determination of the structure of a-Ss at 300 and lOOKc^, where 
the crystalline charge distribution is measured. This study clearly disregard any possible 
structural phase transition and shows that the simple molecule model, in spite of the amount 
of experimental data that it is able to reproduce, should still be improved. 

a) pSs crystals: 

If a-Ss is slowly heated, it shows a phase transition to monoclinic jS-Ss at 369K, which 
melts at 393KEI. The structure has been determined from single crystal dataH at room tem- 
perature; this was possible because, after annealing at 383K, the crystals may be maintained 
at 300K for up to a month. The space group is Cg/i (P2i/c), with six molecules per unit 
cell. Four of the molecules are located at general positions: j:(x,y,z; x,l/2-y,l/2+z) with 
x=0.367, y=0.356 and z=0.072. The other two molecules ( 5 and 6) are located at inversion 
centers: (0,0,0) and (0,1/2,1/2), achieving this symmetry by orientational disorder. Should 
the disorder be static, the averaged symmetry is over all molecules at inversion sites, if the 
disorder is dynamic, the molecules achieve its site symmetry on a time average. In both cases 
the disordered array has a residual entropy of kB/3 ln2 per moleculeH. The experimental 
and lattice parameters of our sample are included in Table I, the molecular orientations are 
as measured in ref.S. 

Fig. 2 includes the calculated configurational energy and volume per molecule, that show 
no discontinuities as a function of temperature. Instead, we find that the reorientation 
of molecules 5 and 6 freezes at low temperatures. This is in accord with heat capacity 
measurements that locate an order- disorder transition transition at 198Ka. Fig. lb shows 
our calculated sample at T=150K, the molecules 5 and 6 are now related by a C2 symmetry 
axes. 

We find that the orientational disorder of molecules at (0,0,0) and (1/2,1/2,1/2), at high 
temperatures, is of the dynamical type. Fig. 3 shows the calculated self correlational func- 



tions for reorientation of the three inertial molecular axes at 370K. The axis perpendicular 
to the molecular plane librates around its equilibrium orientation(z axis in Fig.3), whilst 
the molecule reorientates within the plane. The decay time is about lOps. Furthermore, a 
histogram of molecular orientations, averaged over all molecules and all time steps, shows 
that the initial correlation has been totally lost after 20000 equilibration time steps. For 
T<150K we measure librations of the three axes (with larger amplitudes within the plane), 
but no reorientational diffusion. 

a) 7-S8 crystals: 

7-83 crystals can be grown from solutionallH. The crystal is monoclinic, space group C2/J 
(P2/c), with four molecules per unit celB. This pseudo-hexagonal structure is denser than 
that of a-Sg, at the same temperatureH'H. Our simulations reproduce this structure (Fig. 
Ic), and Fig. 2 shows that at all temperatures this packing has the higher density. Our 
calculated parameters and experimental data, at 300K, are included in Table I. It must be 
noticed that the experimental density, measured at 300K, is 2.19gr/cm^, 5.8% higher than 
that of a-Sgclfi. Nevertheless, in ref.H the X-ray measurement determines a structure with 
volume per molecule higher than that of a-Sg, implying a 2.3% lower density. It could be 
possible that the temperature of the of the X-ray measurement was higher than 300K. 
No other experimental data are available for further comparisons. For this sample we find a 
large discontinuity in lattice parameters for T>400K and molecular diffusion for T>470K. 

d) Phase transitions and the liquid phase: 

There is experimental evidence that the phase transitions are obtained when defects are 
present, if they are not, metastable states can be maintained for daysfil'oH. Our calculations 
do not reproduce the a-(3 transition, but this fact corresponds with the experimental evidence 
that /3-Sg crystals are more easily grown from the melt and not by increasing temperature 
of an a-Sg sample. 

The same evidence exists for the solid-liquid phase transitions. In this case, the exper- 



imental data show that the disorder is generated because Ss molecules start to dissociate 
when the temperature of the sample is near the melting point. 

At even higher temperatures ( 460K) liquid sulfur is known for its liquid-liquid phase tran- 
sition, due to changes in the molecular composition of the elemental sulfur molecules Sn. 

That the disorder is esential to promote the solid-liquid phase transition can be checked 
by simulating a liquid sample of Sg molecules. These simulations should be valid for tem- 
peratures above and near the transition, when the fraction of broken molecules is small. A 
cubic sample of 216 disordered molecules was studied in the range 350-450K, following the 
same procedure as the other samples. In this case the MD algorithm allowed changes in 
the volume of the MD box, but not in its shape. Fig. 4 shows the calculated discontinuity 
in the configurational energy and volume per molecule at 400K, near that measured for a 
(386K) and /? (393K) crystals. For T>400K the sample is liquid, as can be measured by its 
diffusion constant. Comparison of Figs. 2 and 4 show that the solid - liquid transition is 
found if the sample is disordered. 

Conclusions: 

In this paper we study the phase diagram of Ss crystals, at Okbar pressure, as given by a 
simple and flexible molecule model. Our calculations can be also considered a model study 
of cyclo-chain molecules with a torsional potential of a double well type, at difference with 
hydrocarbon chains which show a triple well torsional potential. 

It has to be emphasized that this phase diagram is extremely complex and there are no 
previous calculations on it. The proposed molecular model is extremely simple but, never- 
theless, it gives very good account of many experimental facts. The three crystalline phases 
are stables. The molecular model reproduces the structure, configurational energy and dy- 
namics of a-Ss, /^-Ss and 7-S8 crystals for T>200K. At STP the experimental volume per 
molecule, lattice parameters and molecular array of the three phases compare well with the 
experimental data, including the smallest volumen per molecule for 7-S8 crystals. 
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The calculations also reproduce the orientationally disordered phase of pSs, that we deter- 
mine as djTiamically disordered. At low temperatures we calculate a stable orientationally 
ordered phase. 

The solid-liquid phase transition is found near the experimental value, but only for a disor- 
dered sample. The a- /J-Sg phase transition is not reproduced by either of our two samples. 
Both facts are in accord with the experimental evidence that the transitions are promoted 
by disordeinl, if not, metastable states can be maintained for days. 

Acknowledgement: The authors thank CONICET for the grant PIP 0859/98. 
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Table I: Experimental and calculated lattice parameters: a, b, c (A) and monoclinic 
angle j3 (deg.). Z is the number of molecules in the cell, V (A^) is the volume per molecule. 

See the text for a discussion about the 7-S8 data. 



Phases Z a(A) b(A) c(A) (3 (deg.) Y{k^) 

a-Sg Ref.0 16 10.4646(1) 12.8660(1) 24.4860(3) - 206.046(98) 

(T=300 K) Calc. 10.34(10) 13.20(3) 24.32(7) - 207.3(1.0) 

/3-S8 Ref.i 6 10.778 10.844 10.924 95.8° 211.70 

(T=370 K) Calc. 11.01(11) 10.753(27) 10.79(10) 97.48(07)° 211.30(74) 

7-S8 Ref.i 4 8.442(30) 13.025(10) 9.356(50) 124.98(30)° 210.8(1.3) 

(T=300 K) Calc. 8.03(13) 13.16(14) 8.88(11) 120(3)° 202.5(2.5) 
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Figures: 

Fig. 1: Our calculated crystals: a) orthorhombic a-Sg at 300K, b) monoclinic /S-Sg at 
150K, c) pseudo-hexagonal 7-S8 at 300K. 

Fig. 2: Calculated configurational energy and volume per molecule as a function of 
temperature for all calculated samples. The lines are a guide to the eyes. 

Fig. 3: /S-Sg crystals at 370K: Reorientational self-correlation functions C20 for the 
inertial axes of the two molecules at inversion sites. 

Fig. 4: Phase diagram of a disordered cubic sample, with the same scale of Fig. 2. For 
T>400K, the sample is liquid. 
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Sg is the most stable compound of elemental sulfur in solid and liquid phases, at 
ambient pressure and below 400K. Three crystalline phases of Ss have been clearly 
identified in this range of thermodynamic parameters, although no calculation of its 
phase diagram has been performed yet. a- and 7-S8 are orientationally ordered crys- 
tals while (3-Ss, is measured as orientationally disordered. In this paper we analyze 
the phase diagram of Sg crystals, as given by a simple and flexible molecule model, 
via a series of molecular dynamics (MD) simulations. The calculations are performed 
in the constant pressure- constant temperature ensemble, using an algorithm that is 
able to reproduce structural phase transitions. 

Introduction: 

Elemental sulfur shovv^s a rich variety of molecular forms and therefore of chemical and 
physical properties in the crystalline and liquid phases^'^. Recently, a high pressure super- 
conductor transition has been found, with Tc=17K (at 162GPa), the highest known for any 
elemental solid^. Natural sulfur is of relevance in geophysics, astrophysics, material sciences 
and also of massive industrial use^. Nevertheless, the complex phase diagram and properties 
of its condensed phases are far from well known^'^. 

The most stable sulfur allotrope at ambient temperature and pressure (STP) is Sg, a 
crow-shaped cyclic Sg molecule, stable in solid, liquid and gas phases. 



Orthorhombic a-Ss is the crystalline form stable at STP. Figs la shows its "crankshaft" 
structure with four molecules in the primitive unit cell. The space group is Dl\ (Fddd), with 
16 molecules in this non-primitive orthorhombic cell. If a-Ss is slowly heated, it shows a 
phase transition to monoclinic /J-Ss at 369K (Fig. lb), which melts at 393K^ Nevertheless, 
most tt-Sg crystals do no easily convert to /5-S8, they melt, instead, at 385. 8K^'^. /J-Sg is a 
monoclinic crystal, with six molecules in the primitive unit cell, two of them orientationally 
disordered'''^. A third crystalline allotrope has been observed: 7-S8 (Fig. Ic), that can be 
obtained from solutions of Sg or from its melt^'^. The space group is C2/J (P2/c), with four 
molecules in the pseudo-hexagonal closed-packed unit cell^. The density of this allotrope, 
at STP, is 5.8 % higher than that of a-Sg. 

In spite of its relevance, calculations on the crystalline phase diagram of elemental sulfur 
molecules are scarce. Previous molecular dynamics (MD) simulations of crystalline Ss were 
limited to constant-volume simulations at normal pressure and temperatures of a-Sg^, with 
the only aim to study the intramolecular frequencies obtained by Raman and ir. measure- 
ments. The same scope was followed by numerous works that measured crystalline frequen- 
cies and performed calculations based on lattice dynamics sums in different aproximations^. 
No attempt to reproduce the phase diagram of this molecule, including determination of 
crystalline structures, lattice dynamics, or thermodynamic properties of /5- and 7-83 as a 
function of temperature was intented. Here we study the phase diagram of Ss crystals, as 
given by a simple model of flexible molecule, by performing a series of constant pressure - 
constant temperature MD simulations, at several temperatures and zero pressure. 

The inter- and intramolecular potential model: 

The main problem in this type of calculation is the inter- and intramolecular potential 
model of this flexible molecule. The above mentioned constant volume MD simulations 
of tt-Sg were performed with an intermolecular potential model of the atom-atom type, to 
study the splitting of intramolecular frequencies in a sample of 64 molecules of a -Sg^. In our 
constant temperature - constant pressure MD simulations, these atom-atom parameters are 



slightly changed to improve the fit of our results to the experimental unit cell volume and 
the estimated configurational energy of a-Sg at 300K and Okbar pressure. The parameters 
of our LJ potential model, for S-S interactions are: e= 1.70kJ/mol and (t=3.39 A. 

In a similar way, our intramolecular potential parameters for bending and torsional angles 
of the fiexible Sg molecule are based on those of ref.^. Our molecule model takes into account 
all intramolecular modes that mix with lattice modes and can therefore be relevant for 
structural phase transitions. As calculated in ref.^ and also measured in all crystalline phases 
at all temperatures, three main regions are clearly separated in the vibrational spectra: 
between 10 and ^lOOcm^^ are found the lattice and torsional modes, between 150 and 
250cm~^ the bending modes and from 400 to 500cm~^ the stretching modes, well above in 
energies than the rest of vibrational modes. The bond length is therefore held constant at 
2.O6OA in our calculations. Our fiexible molecule model allows, instead, the mixing of the 
bending and torsional intramolecular motions with lattice modes in a straightforward way. 
The intramolecular potential for bending angles S-S-S is taken armonic 

with a force constant of C/3=25725kB/rad^ and j3q= lOSdeg., k^ is the Boltzman constant. 
The intramolecular potential for torsion (r) angles is a double well: 

V{t) = Ar + BrCos{T) + C^cos^(r) + Dt-cos^{t), 

with ^^=57.192 k^, 5^=738.415 k^, C^=2297.880 k^ and i:'^=557.255 k^. These param- 
eters describe a double well with minima at r=l 98.8deg., and a barrier height of about 
9kJ/mol at r=180deg. The barrier height at r=Odeg. is of 30kJ/mol, out of the range of 
energies explored in these simulations. 

Calculations: 

The phase diagram and dynamical properties of Ss crystals, as given by this molecular 
model, are studied in the (N,P,T) ensemble, at zero pressure, by a series of classical MD 



simulations. The external pressure is considered isotropic and the MD algorithm allows 
volume and shape fluctuations of the sample, in order to balance the internal stresses with 
an external isotropic pressure^. The algorithm controls pressure via an extended system 
which includes as extra variables the MD box parameters. The temperature control of the 
sample follows the approach of Nose^°^^^, which also includes an external variable. The 
equations of motion of these flexible molecules are integrated using the Verlet algorithm for 
the atomic displacements and the Shake algorithm for the constant bond length constraints 
on each molecule^^'^''. The flnal MD algorithm is totally similar to that used in a study of 
black Newton fllms^^. 

As starting point of our simulations, we considered the available structural data at STP. 
The molecules were located at the experimental sites, with the measured orientation but 
with their isolated D4d symmetry. In the course of the simulation the molecules distort to 
attain, on a time average, their site symmetry. Following this procedure, the only problem 
found was with the /3-Ss sample. In this case the calculations were initiated with the LJ 
parameter s multiplied by a factor that in a few (about 100) time steps was varied from 0.1 
to 1. This procedure allowed the molecule to attain its initial distortion in the volume of 
the experimental unit cell. 

In a flrst step we stabilized our samples, at the temperatures quoted in Table I, with a 
constant volume algorithm for 30000 time steps (each one of O.Olps.). Later on, the runs 
in the (N,P,T) ensemble were performed by lowering and increasing the temperature (up 
to 500K) in steps of 25 or 50K. At each point of the phase diagram the samples are sta- 
bilized for 20000 to 30000 time steps and measured in the following lOps. Near the phase 
transitions the stabilization times were increased several times. The sample of a-Sg crystals 
consisted of 3*3*2 orthorhombic cells, a sample of 288 molecules. Some points of the phase 
diagram were recalculated with a sample of 4*4*2 cells (512 molecules). The sample of /J-Ss 
crystals consisted of 4*4*4 cells (384 molecules) and that of 7-S8 crystals of 5*3*5 cells (300 
molecules). 



Results: 

a) a-Sg crystals: 

a-Ss crystallizes, at STP, in the orthorhombic space group F^ (D^/^) with 16 molecules in 
a face-centered unit cell, the primitive cell contains 4 molecules at C2 sites. Fig. la includes 
part of our sample at 300K and shows the structure in the [110] direction. The molecular 
rings lie parallel to the crystallographic c axis and are alternately oriented parallel to (110) 
and (110) crystalline planes. Table I shows the experimental lattice parameters and our 
calculated values, the largest difference is found in the value of b axis, 2 % larger than its 
experimental value. The molecular array and their corresponding orientations follows the 
experimental data. Our calculated configurational energy is -103.5(2)kJ/mol at 300K, from 
which the heat of sublimation can be estimated as lOl.OkJ/mol and compares well with a 
measured value of 101.8(2)kJ/mol^^. Also our calculated pair distribution function g2(r), 
not included here, compares well with neutron measurements of crystalline a-Sg at 300K^^. 

Fig. 2 shows, for ct— , j3— and 7 — S^, the calculated configurational energy and volume 
per molecule as a function of temperature. Our simulations reproduce the experimental 
data of the three crystals at STP. Unfortunately, our a-Sg sample of 288 molecules turns 
out to be unstable for T<200K, and distorts to a monoclinic cell, with molecular locations 
and orientations closely resembling that of Ck-Sg. This result was later on confirmed with a 
larger sample of 512 molecules (4*4*2 orthorhombic cells). Below 200K this large sample 
do not distort to monoclinic, but nevertheless its configurational energy is higher, and its 
volume per molecule is larger, than those calculated for the sample of 288 molecules. Also 
the fiuctuations of the energy and volume, of the sample of 512 molecules, are larger for 
T<200 than their value at 200K. 

The coincidence of the calculated distortion at 200K and the experimental fact that 
crystals of pure Sg are greenish- yellow at 293K and turn white below 193K, made us review 
the experimental data on this phase. Usually the reference quoted for lattice parameters as 
a function of temperature is ref.^^. These values were obtained using X-ray powder diffrac- 
trometry, on different samples and using only three lines. Our calculated powder diffraction 



patterns, not included here, were similar for the orthorhombic and monoclinic samples and 
the scarce experimental data could be not enough to detect the changes. Nevertheless, there 
is a careful single crystal determination of the structure of a-Ss at 300 and 100K^°, where 
the crystalline charge distribution is measured. This study clearly disregard any possible 
structural phase transition and shows that the simple molecule model, in spite of the amount 
of experimental data that is able to reproduce, should still be improved. 

a) /5-S8 crystals: 

If a-Ss is slowly heated, it shows a phase transition to monoclinic /J-Ss at 369K, which 
melts at 393K^. The structure has been determined from single crystal data^ at room tem- 
perature; this was possible because, after annealing at 383K, the crystals may be maintained 
at 300K for up to a month. The space group is C^f^ (P2i/c), with six molecules per unit 
cell. Four of the molecules are located at general positions: l(x,y,z; x,l/2-y,l/2+z) with 
x=0.367, y=0.356 and z=0.072. The other two molecules ( 5 and 6) are located at inversion 
centers: (0,0,0) and (0,1/2,1/2), achieving this symmetry by orientational disorder. Should 
the disorder be static, the averaged symmetry is over all molecules at inversion sites, if the 
disorder is dynamic, the molecules achieve its site symmetry on a time average. In both cas- 
es the disordered array has a residual entropy of 1:^/3 ln2 per molecule^. The experimental 
and lattice parameters of our sample are included in Table I, the molecular orientations are 
as measured in ref.^. 

Fig. 2 includes the calculated configurational energy and volume per molecule, that show 
no discontinuities as a function of temperature. Instead, we find that the reorientation 
of molecules 5 and 6 freezes at low temperatures. This is in accord with heat capacity 
measurements that locate an order-disorder transition transition at 198K''. Fig. lb shows 
our calculated sample at T=150K, the molecules 5 and 6 are now related by a C2 symmetry 
axes. 

We find that the orientational disorder of molecules at (0,0,0) and (1/2,1/2,1/2), at high 
temperatures, is of the dynamical type. Fig. 3 shows the calculated self correlational func- 
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tions for reorientation of the three inertial molecular axes at 370K. The axis perpendicular 
to the molecular plane librates around its equilibrium orientation(z axis in Fig. 3), whilst 
the molecule reorientates within the plane. The decay time is about lOps. Furthermore, a 
histogram of molecular orientations, averaged over all molecules and all time steps, shows 
that the initial correlation has been totally lost after 20000 equilibration time steps. For 
T<150K we measure librations of the three axes (with larger amplitudes within the plane), 
but no reorientational diffusion. 

a) 7-S8 crystals: 

7-S8 crystals can be grown from solutions^'^. The crystal is monoclinic, space group C2/J 
(P2/c), with four molecules per unit cell^. This pseudo-hexagonal structure is denser than 
that of ct-Ss, at the same temperature^'^. Our simulations reproduce this structure (Fig. 
Ic), and Fig. 2 shows that at all temperatures this packing has the higher density. Our 
calculated parameters and experimental data, at 300K, are included in Table I. It must be 
noticed that the experimental density, measured at 300K, is 2.19gr/cm^, 5.8% higher than 
that of tt-Sg^'^. Nevertheless, in ref.^ the X-ray measurement determines a structure with 
volume per molecule higher than that of ct-Ss, implying a 2.3% lower density. It could be 
possible that the temperature of the of the X-ray measurement was higher than 300K. 
No other experimental data is available for further comparisons. For this sample we find a 
large discontinuity in lattice parameters for T>400K and molecular diffusion for T>470K. 

d) Phase transitions and the Hquid phase: 

There is experimental evidence that the phase transitions are obtained when defects are 
present, if they are not, metastable states can be maintained for days^'^'^. Our calculations 
do not reproduce the a-/? transition, but this fact corresponds with the experimental evidence 
that /J-Sg crystals are more easily grown from the melt and not by increasing temperature 
of an ct-Ss sample. 

The same evidence exists for the solid-liquid phase transitions. In this case, the exper- 



imental data shows that the disorder is generated because Ss molecules start to disociate 
when the temperature of the sample is near the melting point. 

At even higher temperatures ( 460K) liquid sulfur is known for its liquid-liquid phase tran- 
sition, due to changes in the molecular composition of the elemental sulfur molecules Sn. 

That the disorder is esential to promote the solid-liquid phase transition can be checked 
by simulating a liquid sample of Sg molecules. These simulations should be valid for tem- 
peratures above and near the transition, when the fraction of broken molecules is small. A 
cubic sample of 216 disordered molecules was studied in the range 350-450K, following the 
same procedure as the other samples. In this case the MD algorithm allowed changes in 
the volume of the MD box, but not in its shape. Fig. 4 shows the calculated discontinuity 
in the configurational energy and volume per molecule at 400K, near that measured for a 
(386K) and /? (393K) crystals. For T>400K the sample is liquid, as can be measured by its 
diffusion constant. Comparison of Figs. 2 and 4 show that the solid - liquid transition is 
found if the sample is disordered. 

Conclusions: 

In this paper we study the phase diagram of Sg crystals, at Okbar pressure, as given by a 
simple and flexible molecule model. Our calculations can be also considered a model study 
of cyclo-chain molecules with a torsional potential of a double well type, at difference with 
hydrocarbon chains which show a triple well torsional potential. 

It has to be emphasized that this phase diagram is extremely complex and there are no 
previous calculations on it. The proposed molecular model is extremely simple but, never- 
theless, it gives very good account of many experimental facts. The three crystalline phases 
are stables. The molecular model reproduces the structure, configurational energy and dy- 
namics of tt-Sg, /3-Sg and 7-Sg crystals for T>200K. At STP the experimental volume per 
molecule, lattice parameters and molecular array of the three phases compare well with the 
experimental data, including the smallest volumen per molecule for 7-Sg crystals. 
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The calculations also reproduce the orientationally disordered phase of /J-Ss, that we deter- 
mine as dynamically disordered. At low temperatures we calculate a stable orientationally 
ordered phase. 

The solid-liquid phase transition is found near the experimental value, but only for a disor- 
dered sample. The a- /J-Sg phase transition is not reproduced by any of our two samples. 
Both facts are in accord with the experimental evidence that the transitions are promoted 
by disorder^, if not, metastable states can be mantained for days. 
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Table I: Experimental and calculated lattice parameters: a, b, c (A) and monoclinic 
angle /? (deg.). Z is the number of molecules in the cell, V (A^) is the volume per molecule. 
See the text for a discussion about the 7-S8 data. 



Phases Z a(A) b(A) c(A) /3 (deg.) V(A^) 

a-Sg Rei}^ 16 10.4646(1) 12.8660(1) 24.4860(3) - 206.046(98) 

(T=300 K) Calc. 10.34(10) 13.20(3) 24.32(7) - 207.3(1.0) 

P-Sg Ref.2 6 10.778 10.844 10.924 95.8° 211.70 

(T=370 K) Calc. 11.01(11) 10.753(27) 10.79(10) 97.48(07)° 211.30(74) 

7-S8 Ref.6 4 8.442(30) 13.025(10) 9.356(50) 124.98(30)° 210.8(1.3) 

(T=300 K) Calc. 8.03(13) 13.16(14) 8.88(11) 120(3)° 202.5(2.5) 
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Figures: 

Fig. 1: Our calculated crystals: a) orthorhombic a-Sg at 300K, b) monoclinic /J-Sg at 
150K, c) pseudo-hexagonal 7-S8 at 300K. 

Fig. 2: Calculated configurational energy and volume per molecule as a function of 
temperature for all calculated samples. The lines are a guide to the eyes. 

Fig. 3: /J-Ss crystals at 370K: Reorientational self-correlation functions C20 for the 
inertial axes of the two molecules at inversion sites. 

Fig. 4: Phase diagram of a disordered cubic sample, with the same scale of Fig. 2. For 
T>400K, the sample is liquid. 
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